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ABSTRACT 

We present results of the near-infrared (IR) spectroscopy of six quasars whose 
redshifts range from 0.158 to 1.084. Combined with the satellite ultraviolet data, 
the relative line strengths of O I A1304, O I A8446, O I A11287, and the near-IR 
Ca II triplet are given. In addition, the corresponding O I line strengths measured 
in normal Seyfert Is and narrow-line Seyfert Is are collected from the literature. 
These lines are thought to emerge from the same gas as do the Fe II lines, 
so they are good tracers of the Fe II emission region within a broad emission 
line region (BELR) in active galactic nuclei (AGNs). In order to reveal the 
physical condition within the relevant emission region, we performed photoionized 
model calculations and compared them to the observations. It suggests that a 
rather dense gas with density nn ~ 10 115 cm" 3 is present at an outer portion of 
the BELR, illuminated by the ionizing radiation corresponding to an ionization 
parameter U ~ 10~ 25 and is primarily responsible for the observed O I, Ca II, 
and Fe II lines, based on the resemblance of their profiles. The three O I lines 
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are proven to be formed through Ly/3 fluorescence and collisional excitation. We 
also show that the A1304 bump typically observed in AGN spectra consists of 
the comparable contributions of O I and Si II multiplets, and we discuss the 
origin of such a strong Si II emission. The results are interpreted in the context 
of the locally optimally emitting cloud (LOC) scenario to find the plausible gas 
distribution within the BELR as a function of distance from the central source 
and density. 

Subject headings: galaxies: active — galaxies: individual (3C 273, QSO B0850+440, 
3C 232, QSO J1139-1350, PG 1148+549, PG 1718+481) — galaxies: quasars: 
emission lines — galaxies: Seyfert — line: formation 



INTRODUCTION 



The emission lines from singly ionized iron are one of the most prominent features in 
ultraviolet (UV) to optical spectra of many AGNs. They are considered to arise from the 
partly ionized gas present in a BELR. They dominate the heating and cooling of the emission 
region through a huge number of electronic transitions characteristic of the Fe + atom, which 
results in the Fe II "pseudocontinuum" observed in the AGN spectra. 

It is important to understand the physics of the Fe II line formation for several rea- 
sons. First, since the Fe II lines carry a large amount of energy away from the BELR, 
it provides useful information about the energy budget of the emitting gas. Second, mea- 
surements of the Fe abundance in high-redshift quasars through thei r Fe II emission as 
a function of cosmic tim e could verify some cosmological parameters (IHamann fe Ferland 
19931 : lYoshii et al.l Il998l ). Many observations have been devoted to t he Fe II flux mea- 



surements in hiRh-r e dshift quasars for th i s purpose in recent years (e.g.. lElston et al. 



Kawara et al. 



2003: 



1996 



Maiolino et al 



Verner et al 



2003 



Dietrich et al.l 12002 . 12003 : llwamuro et all I2002L 12004 : iFreudling et al. 



1994 



2003 ). However, as has been pointe d out by several authors (e.g. 



Baldwin et al.l l2004t iTsuzuki et al.l 120071 ) , the Fe II line strengths could 



be strongly affected by some physical, nonabundance factors of the emitting gas, such as 
density and incident-ionizing radiation flux, so that a priori knowledge of them is definitely 
needed in order to deduce the Fe abundance. 



There have b een many studies of Fe II emission in AGNs. iNetzer fe Willa (119831 ) and 
Wills et al.l (119851 ) modeled the Fe II emission of quasars using a model Fe + atom with energy 
levels up to ~10 eV and typical BELR p arameters in order to compare the m with obser 



vations. Their study was followed by, e.g., ISigut fe Pradhanl (119981 . 120031 ) and IVerner et al. 
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( 119991 . 120031 ) , who made significant imp r oveme nts both on the atomic data and on the phys- 
ical processes included. iBaldwin et al.l (120041 ) calculated the Fe II and other emission-line 
strengths in a wide range of gas physical par ameters and integrat ed them to reproduce the 
observed spectra in the context of the LOC (IBaldwin et al.lll995l ) scenario. However, these 
studies using the observed shape and/or strength of the Fe II pseudocontinuum have several 
serious problems. From an observational viewpoint, it is extremely difficult to isolate the 
Fe II pseudocontinuum from the UV to optical spectrum of AGNs, which is usually made up 
of several continuum component s besides emiss i on lin es, such as a power-law continuum and 
a Balmer continuum. Recently, iTsuzuki et al.l (120061 ) presented quasar spectra with a wide 
wavelength coverage, extending from 1000 to 7300 A, thus providing a detailed discussion 
on this issue. From a theoretical viewpoint, it is still not confirmed whether the included 
number of energy levels of the Fe + atom and excitation mechanisms are sufficient to repro- 
duce the physical processes occurring in AGNs. Specifically, it is widely known that none 
of th e current photoionize d models can reproduce the observed strong optical Fe II emitters 



(e.g., iCollin fe Jolyl |2000| ) while the recent reverberation mapping results indicate a direct 
connect ion between the optical Fe II e mission and the continuum radiation from the central 
source (jVestergaard fe Peterson! 120051 ). 



Given the difficulties raised above, one approach to investigating the physical condition 
within the Fe II emission region is to use the lines emitted by simple atoms in the same region. 
Rodriguez- Ardila et al.l (l2002al ) found that near-IR O I, Ca II, and Fe II lines have very 
similar profiles in the narrow-line Seyfert Is (NLSls), indicating that these low-ionization 
lines are arising from the same portion of the BELR. They also found a good correlation 
between the near-IR and optical Fe II emission, which implies a direct link between them 
in the sense that they are produced by the same excitation mechanisms. The correlation 
between O I and Ca II line widths and betw een Ca II and Fe II line strengths has also 
been reported for Seyfert Is (Syls) and NLSls (|Perssonlll988l ). Note that they are a natural 
consequence of similar ionization potentials of the relevant ions, i.e., 16.2 eV for Fe II, 13.6 
eV for O I, and 11.9 eV for Ca II. Thus, O I and Ca II lines are good tracers of the interesting 
and important Fe II line-emission region in AGNs. 

The first ext ensive study of the physical condition within an O I emission region in AGNs 
was presented by iGrandil (Il980l ). He found that O I A8446 lacks the narrow component that 
characterizes other permitted lines, leading to the c onclusion that the li ne is pur e ly a B ELR 
phenomenon. Such a conclusion is also reported in iRiffel et al.l (120061 ) . IGrandil (Il980l ) also 
suggested that O I A8446 is produced by Ly/3 fl uorescence, which w a s later confirmed by the 



observa tion of I Zw 1, the prototype NLS1, by iRudy et al.l (119891 ). iRodriguez- Ardila et al. 



(l2002bl ) compiled the UV and near-IR O I line strengths of NLSls (including I Zw 1) and of 
Syls, to investigate their flux ratios, and found that there must be an additional excitation 
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mechanism of I A8446 — - b esides Ly/3 fluorescence — which they concluded is collisional 
excitation. IMatsuoka et al.l (120051 ) found that the collisional processes are at work to a 
similar extent in quasars and suggested that the gas density in the relevant line emission 
region is not signif i cantly different among NLSls, Syls, and quasars. On the other hand, 



Ferland fe Perssonl (119891 ) computed the strengths of various emission lines, including the 
UV and near-IR Ca II lines, and found that a very large column density (~ 10 25 cm -2 ) 
is needed in order to reproduce the observed relative strengths of these lines. However, 
previous studies have not yet revealed some important physical parameters of the emitting 
gas quantitatively, mainly due to the lack of sufficient data or precise modeling. 

Here we present a new near-IR observation of six quasars with redshifts up to 1.084. 
Combined with the satellite UV data, the relative strengths of O I A1304, O I A8446, O I 
A11287, and the near-IR Ca II triplet (A8498, A8542, and A8662) are given. The purpose 
of the observation is twofold: to give a quantitative insight into the physical condition 
within the partly ionized region where O I, Ca II, and Fe II lines are emitted, by comparing 
the observations with precise model calculations; and to extend the observed O I and Ca II 
samples to higher redshifts than the previous measurements (z < 0.18) and to include quasars, 
which are important in the context of the Fe abundance investigation in high-redshift quasars. 
The observations are presented in SJ2], and compared model calculations are given in §3j A 
discussion of the results and their implications appear in S]D Finally our conclusions are 
summarized in 35l 



2. OBSERVATIONS 
2.1. Target Selection 



Target quasars were selected from a UV spectral atlas of AGNs presented by Evans fe Koratkar 



( 120041 ) . The strongest constraint on the target selection comes from their redshifts; O I A8446, 
O I A11287, and the near-IR Ca II triplet should fall within the atmospheric windows. It 
significantly reduced a number of possible targets, since these lines are too separated to be 
observed in the single near-IR band (J, H, or K). From the quasars satisfying this criterion, 
we selected those having the highest signal-to-noise ratios (S/Ns) around O I A1304 in the 
corresponding UV spectra, which constitute our final target list. The characteristics of the 
observed quasars are summarized in Table [H 
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2.2. Near-Infrared Observation 



A near-IR observation was carried out at the United Kingdom Infrared Telescope 
(UKIRT) using the UIST (UKIRT 1-5 fim Imager Spectrometer; iRamsay Howatt et al. I l2004h 



on 2006 February 24. We adopted the standard UIST observing templates for faint point- 
source spectroscopy and for point-source imaging. An observational journal is given in Table 
[2] for the spectroscopy and in Table [3] for the imaging. The spectroscopic observation in- 
cludes a flat, arc, and rationing star for each target, taken right before the target observation. 
The rationing stars were selected from the F4 - F9 type dwarf bright stars with moderately 
strong H I recombination lines and weak atomic lines, which are well suited for rationing 
low-resolution spectra. They were observed at air masses similar to the corresponding tar- 
gets. The targets and the rationing stars were observed by nodding the telescope in such 
a way that they slid +6" and then —6" along the slit. A 0.6" slit yielded an instrumental 
resolution of ~750 km s _1 , while the pixel scale is 0.12". The imaging observation consists 
of dark, flat, and standard stars taken right before the target observation, while the self-flats 
were used for some faint sources. The standard stars were taken from the UKIRT Faint 



Standards (IHawarden et al.ll200ll ; iLegett et al.ll2006l ) so that they could be observed at air 



masses similar to the targets. The targets and the standard stars were observed at five po- 
sitions on the detector, offset by 20" for the targets and by 10" for the standard stars. The 
sky condition was photometric throughout the night, and the average seeing was 0.6". 

We used the standard method to reduce the data, including removal of dark, flat- 
fielding, and sky subtraction using the reduction software IRAFj^J Wavelength calibration 
was achieved using argon arc-lamp spectra. In order to calibrate the sensitivity variation 
along the dispersion axis, the rationing-star spectra — from which the intrinsic H I recombi- 
nation absorption lines had been removed — were compared with the blackbody spectra with 
corresponding temperatures, i.e., 6500 K for F4 V type stars, 6250 K for F7 V, and 6000 
K for F8 V and F9 V. Finally, photometric calibration was obtained by using the imaging 
data. Measured J- or if -band magnitudes of the quasars are listed in the last two columns 
of Table [TJ Statistical errors of the photometry were less than 0.01 mag for all targets. 

We show the reduced spectra along with the sensitivity curves in Figures [TJ - [2j The 
redshifts measur ed in the near-IR spectr a were in good agreement with those measured in 



the UV spectra (jEvans fc Koratkarl 120041 ) 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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2.3. Ultraviolet Data 



Ultraviolet spectra of the quasars were obtained from 1991 January to 1993 July using 
the Faint Object Spectrograph (FOS) o nboard the Hubble Space Telescop^ (HST). In the 
current work, we used those reduced by lEvans &: Koratkarl (12004 ) . They extensively recal- 
ibrated all pre-COSTAR (Corrective Optics Space Telescope Axial Replacement) archival 
FOS UV and optical spectrophotometry of AGNs with the latest algorithms and calibration 
data. If multiple observations of the same source were available, they were combined after 
being scaled to each other in such a way that the resultant spectrum would have the highest 
possible S/N. The UV spectra around O I A1304 used in this work are shown in Figure [31 



2.4. Emission-Line Measurements 



Emission lines in the UV and the near-IR spectra were manually identified, and their 
intensities were measured by summing up all the flux above the local continuum, which 
was placed by the interpolation method using the observed fluxes in the continuum win- 
dows at both sides of the lines. The continuum windows were carefully selected to not be 
co ntaminated by wea k emi ssion lines based on the high- S /N spectra of NLSls presented 
by lLaor et al.l (ll997bl ) and iRodriguez-Ardila et al.l (j2002al ) . A Voigt profile was fitted to 
continuum-subtracted spectra by the least y 2 metho d to measure line widths for some lines 
(see below). As described by iGlikman et al.l (120061 ). a single-component Voigt profile is a 
better choice than a mult i- component Gaussian profile when fitting lines with a low S/N. 
This is the case for some of the measured lines in our spectra. When several lines were 
blended, we fitted a mult i- component Voigt profile to the continuum-subtracted spectra in 
order to deblend the feature. Then the total flux was divided into each component, according 
to the fitted profiles. 

Although the A1304 bump results from a blending of an O I triplet at A1302, A1305, 
and A1306, and a Si II doublet at A1304 and A1309, deblending of these components is 
hopeless, considering the broad line widths of quasars. Thus, we regarded this feature as a 
single component in the flux-measurement procedure. Its constitution is closely discussed 
in P~21 On the other hand, when fitting blended O I A8446 and the Ca II triplet (A8498, 
A8542, and A8662), we fixed their relative wa yelength centroid s and relative intensities of 
the Ca II triplet to the values found in I Zw 1 (IRudy et al.ll2000l ) in order to avoid too many 



2 This publication is based on observations made with the NASA/ESA Hubble Space Telescope, obtained 
from the data archive at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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parameters being used in the fitting. This procedure was proven to work very well, since 
the blueward part of O I A8446 is almost free from contamination of the Ca II components. 
In some objects, there is an unidentified weak emission line at the blueward part of the 
blend feature, or the redward part falls in the atmospheric absorption bands, which limits 
the available wavelength range used in the fitting and flux measurements. In these cases, 
the correction was applied to the flux falling in outside the flux-measured wavelength range 
according to the fitted Voigt profiles. The fits for the relevant O I and Ca II lines in 3C 273 
are shown in Figure HI 

The measured line fluxes are summarized in Table HI The quoted errors were estimated 
by the A% 2 test, which reflect the uncertainties in determining the underlying continuum 
level (i.e., the S/N in the continuum windows) and in deblending the blended feature. 



2.5. Variability and Reddening Effects 



It is widely known that quasars are highly variable. Since the dates of the UV obser- 
vations are separated from that of the near-IR observation by more than 15 yr, a possible 
variability effect should be carefully examined before comparing the UV and the near-IR 
spectra. In fact, the observed quasars show brightness variations in the J or H band by as 
large as 0.6 mag from the dates of the Two Micron All Sky Surve}]^ (2MASS) observations as 
shown in Table [TJ However, amounts of the emission-line variabilities corresponding to these 
continuum variations are rather uncertain. Monitoring observations of the well-known vari- 
able Seyfert NGC 5548 revealed that when the UV continuum varied its brightness by about 
300%, the flux v ariation in Mg II, one of the l ow-ionization lines su ch as O I and Ca II, was 
only a bout 10% ([Dietrich fe Kollatschnylll995l ). On the other hand. IVestergaard fe Peterson 
( 120051 ) showed that the variability amplitude of optical Fe II emission is 50-75% of that 
of H/3, whose variation is comparable to those in the continuum, and UV Fe II emission 
varies with an even larger amplitude. The response of the individual lines to the continuum 
variation depends o n their formation processes and gas distribution within the BELR (e.g., 
O'Brien et al.lll995l ). whose quantitative estimates are far beyond the scope of this paper. 
We give notes as to how the possible variability could affect our conclusions in the following 
sections, when needed. 



Another source of errors in comparing the UV and the near-IR spectra is the extinc- 



3 This publication makes use of data products from 2MASS, which is a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded 
by the National Aeronautics and Space Administration and the National Science Foundation. 
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tion by dust both in the Gal axy and in the quasar host galaxies. We adopted the Galactic 
extinction values presented by lSchlegel et al.l (119981 ). which are derived from the far-IR obser- 
vations in two satellite missions, the Infrared Astronomical Satellite (IRAS) and the Cosmic 
Background Explorer (COBE). Their extinction map has a relatively fine resolution of 6.1' 
and an accuracy of 16%. Corresponding d ereddening of the spectra was achieved by using 
the Galactic extinction curve presented by iPeil (119921 ). 



On the other ha n d, m uch less has been confirmed about the internal reddening of 



quasars. iCheng et al.l (1199ll ) found that UV continuum slopes for a wide variety of AGNs 
were essentially the same to within their me asuring errors, corr e spond ing; to a reddening 
amount of Eb-v < 0.1 mag. In the optical, iDe Zotti fc Gaskelll (119851 ) showed that both 
the broad emission lines and continua of radio-quiet Seyfert galaxies typically have a red- 
dening of -Sts-v = 0- 25 mag. Assu ming the Sm all Magellanic Cloud (SMC) extinction curve 
( iHopkins et al.ll2004l ) presented by IPeil (119921 ). this translates to -E1304-8446, the color excess 



between 1304 and 8446 A, of 3.58 mag. Hence, the measured UV to near-IR line ratios 
should be regarded as the lower limits on their intrinsic values, which could be an order 
of magnitude or more larger than observed. As in the case of the variability, the possible 
changes in our conclusions due to the internal reddening effect are noted in the following 
sections. 



3. COMPARISON WITH MODELS 



In this section, we aim to investigate the physical condition within the O I and Ca II 
emission region by comparing the observed line strengths with model predictions. It is 
becoming a common understanding that the BELR is composed of gas with widely dis- 
tributed physical parameters and that each emission line arises from its preferable envi- 
ronment ([Baldwin et al.lll995l ). In order to reveal the distribution functions of gas in such 
a complex system, it is crucial to constrain which region in the parameter space actually 
dominates observed emission of the specific lines. 



3.1. Observational Constraints 

The observed equivalent widths (EWs) of O I A8446 and photon flux ratios of O I A11287 
and the near-IR Ca II triplet to O I A8446 are used to constrain the models. Hereafter, they 
are expressed as EW (O I A8446), O I n(A11287)/n(A8446), and n(Ca ll)/n(0 I A8446), 
respectively. These emission lines have been revealed to arise from the same portion of the 
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BELR among the various lines observed (see JQ), which significantly simplifies the interpre- 
tation of the observed quantities. They are summarized in Table |5] (see the quasar section) 
for the quasars obs e rved in this work and for the quasar PG 1116+215 (z = 0.176) studied 



by iMatsuoka et al.l ( 20051 ) . A representative wavelength of 8579 A was used to convert total 



flux of the Ca II triplet to the photon number flux. 

The observed values of O I n(A11287)/n(A8446) are fairly similar from quasar to quasar 
except for 3C 232. The mean value is 0.46 ± 0.18, where the quoted error reflects a scatter 
among the quasars. It should be an excellent indicator of the emitting gas physics, since 
several important but uncertain parameters could be nearly canceled out, which include 
chemical composition and the covering fraction of the line-emitting gas seen from the central 
continuum source. The reddening effect on the ratio would also be small. That observed in 
3C 232 is significantly different from other quasars, which is apparently due to the extreme 
weakness of O I A8446. It is noted that Ca II is also very weak in 3C 232. Since there are no 
distinct features in the sensitivity curve around the relevant wavelength range, the weakness 
is not a consequence of the instrumental or the atmospheric absorption effects. On the other 
hand, no appropriate models considered below could reproduce such a large value of the O I 
n(A11287)/n(A8446) ratio. Hence, we concluded that some unusual effects concealed a large 
fraction of emission lines around 8500 A and decided to exclude this quasar from the sample 
to be compared with model calculations below. 

EW (O I A8446) represents the efficiency of reprocessing the incident continuum radia- 
tion into the line emission. It also depends on the gas covering factor and abundance of the 
relevant element (oxygen in this case). Models should reproduce the observed values, EW 
(O I A8446) > 10 A, with reasonable values of these parameters. 

The Ca II triplet is observed to have a strength comparable to O I A8446. The ratio 
between them, predicted in model calculations, also depends on the assum ed chemical com- 



position (Ca/O in this case) of the emitting gas. lHamann Sz Ferlandl (119991 ) showed that the 
Mg/O abundance ratio varies from the solar value within a factor of 2 through the course of 
galaxy evolution in either the solar neighborhood or the giant elliptical model. Since Ca and 
Mg belong to the a-elements so that their production rates should not significantly differ 
from each other, the Ca/O ratio could also deviate from the solar value by a similar amount. 



3.2. Model Calculations 



A partial Grotrian diagram of an O I atom is given in Figure [5j In the simplest case, O I 
A11287, A8446, and A1304 are formed through Bowen resonance-fluorescence by Ly/5, or Ly/3 
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fluorescence; Ly/3 photons pump up the ground-term (2p 3 P) electrons to the excited 3d 3 D° 
level, from which they cascade down to the ground term through 3p 3 P and 3s 3 S°, producing 
the A11287, A8446, and A1304 line photons. Note that the pumped 2p 3 P - 3d 3 D° transition 
energy corresponds to 1025.77 A so that it falls within a Doppler core of Lyp (1025.72 A) 
for gas at a temperature of 10 4 K. As a result, the photon flux ratios between these lines 
should be unity. However, there could be other excitation mechanisms that alter the ratios, 
which are collisional excitation, continuum fluorescence, and recombination. iKwan fc Krolik 
(119811 ) reported that trapping of A1304 becomes important in sustaining the population at 
3s 3 S° level with the large column density typical of the BELR cloud, leading to collisional 
excitation of the A8446 emission. If A8446 also became optically thick, an additional A11287 
excitation would follow. Continuum fluorescence of the ground-term (2p 3 P) electrons could 
also produce these lines, as well as other lines, such as A13165 (see Fig. [5]). This mechanism 
is know n to be an imp or tant contributor t o the A8446 emission observed in some Galactic 
objects (Grandi 19 75 al ibi; iRudy et~ab 1991 ). While recombination is another source of the 
O I line emissions, its contribu tion is known to be negligible in the BELR gas (lGrandilll980l ; 



Rodriguez- Ardila et al.ll2002bf ). 



We performed model calculations in the fra mework of the phot oionized BELR gas using 
the photoionization code Cloudy, version 06.02 (IFerland et al.lll998l ). The chosen form of the 
incident continuum shape is a combination of a f u oc u auv exp (—hv/kT cnt ) UV bump with an 
X-ray power law of the form /„ oc v a * spanning 13.6 eV to 100 keV. The UV bump is also cut 
off in the infrared with a temperature kTiR = 0.01 Ryd, corresponding to 9.1 /im. The UV 
and X-ray continuum components are combined using a UV to X-ray logarithmic spectral 
slope a ox , which is defined by f v {2 kev)// i/ (2500 A) = 403.3 Qox . We adopted a parameter 
set of [ Tent, q 11v , ay, aw = ( 1 .5 x 10 5 ) K, —0.2, —1.8, —1.4]. It is essentially the same as 
those of iTsuzuki et al.l (120061 . 120071 ) who deri yed these parameter values by combining the 



results of the thermal accretion disk models (IKrolik fc Kallmanl Il988l ; iBinette et al.l 11989 



Zheng et al.lll995l ) and the observed spectrum of the quasar PG 1626+554 which has typical 



spectral properties among their 14 quasar samples. 

The BELR gas was modeled to have a constant hydrogen density ?t.h (cm -3 ) and exposed 
to the ionizing continuum radiation with a photon flux of $ (s -1 cm -2 ). The incident-ionizing 
continuum flux is expressed with an ionization parameter in the calculation U = $/(nn c), 
where c is the speed of light, since predicted line strengths tend to be similar along the same 
values of U. We performed the calculations with the U) sets in a range of 10 7 cm -3 
< % < 10 140 cm -3 and 10 -50 < U < 10 00 stepped by 0.5 dex, which should more than 
cover the parameter space of the BELR gas. We began with a BELR gas with = 10 23 
cm -2 and f tur b = km s -1 , where iVR- is the gas column density and v tU rb the microturbulent 
velocity. Chemical composition was assumed to be solar. We hereafter refer to this model 
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as the standard model, or model 1; input parameters are summarized in Table [61 



In the calculations, the O I atom was treated as a 6 electron energ y-level system includ- 
ing 3s 3 S°, 4s 3 S°, 2p 3 P, 3p 3 P, 4p 3 P, and 3d 3 D° levels (see Fig. EJ iFerlandlliooi l. Note 
that the contribution from recombination is not included in the predicted line strengths, 
since the effective recombination coefficients for the O I lines are not known for densities 
as large as those in the BELR (G. J. Ferland 2006, private communication). Howe ver, it is 
expec ted t o be negligible compared to oth er excitation mechanisms as reported by iGrandi 
(Il980l ) and lRodriguez-Ardila et al.l (j2002bl ). They showed that the observed strengths of O I 
A7774, the quintet counterpart to O I A8446, were less than 10% of those of A8446 (see Table 
IHJ), while the ratio between them should be A7774/ A8446 ~_l-7 i n the pure recombination 



case based solely on the relative statistical weights (IGrandil Il980l ) . On the other hand, we 
found that the ratio should be A7774/A8446 ~ 7.2 in the pure recombination case using 
the effective recombination coefficients in the low-densi t y lim it (appropriate for the normal 
Galactic gaseous nebulae) presented by lPequignot et al.l (1199 ll ). Even if we assumed that the 
observed A7774 was exclusively produced by recombination, its contribution to the A8446 is 
less than a few percent. 



3.3. Results 

3.3.1. The Standard Model 

We show the calculated O I photon flux ratios n(A11287) /n(A8446) and n(A1304) /n(A8446) 
in Figure [6] as a function of the gas density % and the ionization parameter U. Note that 
the incident-ionizing continuum flux $ increases toward the direction shown by the arrows 
in the figure. Both ratios are almost unity at the canonical BELR parameters (tih, U) = 
(10 10 cm " 3 , 10~ 2 ), which refer to the region where high-ionization lines (such as C IV) are 



emitted ([Davidson fc Netzerlll979l ). indicating that Ly/3 fluorescence is actually the main for- 
mation mechanism of these O I lines. As the density becomes large relative to the canonical 
parameter, both ratios decrease due to the collisional enhancement of the A8446 photons. 
On the other hand, continuum fluorescence becomes important as incident continuum flux 
is significantly decreased, due to a significant suppression in Ly/3 formation. It results in 
A 1304 photons not accompanied by those of A8446, as well as A8446 not accompanied by 
A11287, so that n(A11287) < n(A8446) < n(A1304). The results of our observation, O I 
n(A11287)/n(A8446) = 0.46 ± 0.18, clearly indicate that Ly/3 fluorescence is not the only 
excitation mechanism of these lines. The excess A8446 photons relative to A 11287 should be 
produced through either collisional excitation or continuum fluorescence, depending on the 
physical condition within the emitting gas. 



-12- 



The two excitation mechanisms can be distinguished by their continuum reprocessing 
efficiency, which converts the incident continuum radiation into the relevant line emission, 
through the observed EW (O I A8446). We plot the calculated values of them in Figure 
[7J (top right), assuming a covering factor of 1.0, as well as the O I n(A11287)/n(A8446) 
ratios for reference (top left). Since the covering factor in actual quasars may be much 
smaller, all the parameter sets satisfying EW (O I A8446) > 10 A are consistent with the 
observation. However, it is noted that the covering factor should not be much less than 0.1; 
in such a situation, the predicted EW (O I A8446) falls well below those observed on the 
whole parameter plane. It is also true in other models described below. Figure [7] shows 
that the low density gas in which continuum fluorescence dominates the O I line formation 
predicts O I A8446 at the level of EW (O I A8446) < 1 A, which can hardly reproduce the 
observation. On the other hand, there is a large area consistent with both the observed O I 
n(A11287)/n(A8446) ratios and EW (O I A8446) at the high-density side where collisional 
excitation is at work. Thus, we concluded that the observed O I lines are produced by a 
combination of Ly/3 fluorescence and collisional excitation. 

We show the calculated n(Ca ll)/n(0 I A8446) ratios in Figured (bottom left) to further 
constrain the model parameters. The fact that these two emissions are observed to have 
comparable strengths, combined with the results from the O I line strengths, finely restricts 
the gas density to be around = 10 11 ' 5 cm -3 . However, the model predictions could vary 
within a factor of 2, depending on the assumed chemical composition as described in §3.11 It 
also depends on the assumed column density N-r and increases by up to a factor of 3 when 
it is changed from Nn = 10 23 to 10 25 cm -2 , as shown in Figure [HI Thus, we required that 
the predicted Ca II/O I A8446 ratio should fall in the range 0.1 < n(Ca II) / n(0 I A8446) 
< 5.0. This constraint still effectively excludes the high-density (tih > 10 12,0 cm -3 ) and 
low ionization parameter (U < 10 -2 ' 5 ) models in which the Ca II emission is significantly 
overpredicted. 

The top left, top right, and bottom left panels of Figure UJ are shown in a single plot 
in the bottom right panel of the same figure. The constraints from the observed values of 
O I n(A11287)/n(A8446), EW (O I A8446), and n(Ca II)/n(0 I A8446) are shown by thick 
solid lines, thin solid lines, and the shaded area, respectively. The models consistent with all 
the constraints are marked with diamonds, which are wh = 10 11 ' 5 cm" 3 , and U = 10~ 2,5 and 
10~ 3 ' . However, since the covering factors in actual quasars are known to be ~0.1, we favor 
the models predicting EW (O I A8446) > 100 A despite the formal constraint of > 10 A. 
Thus, we concluded that the best-fit parameter set to the observations is (nn, U) = (10 11,5 
cm" 3 , 10~ 30 ) in the standard model. The predicted line strengths in this specific model are 
summarized in Table UJ 
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Next we investigated the line strengths as a function of gas column density in the range 
Nn = 10 17 — 10 25 cm" 2 . A calculation was performed for the gas with Nn = 10 25 cm -2 , while 
other input parameters are identical to the standard model with (ng, U) = (10 115 cm -3 , 
10- 3 - ), and the cumulative ratios of O I n(A11287)/n(A8446), O I n(A1304)/n(A8446), and 
n(Ca II)/n(0 I A8446), as well as the cumulative flux of O I A8446, are plotted against the 
column density measured from the illuminated surface in Figure [HJ The O I A8446 emission 
is mainly built up before the column density Nr = 10 23 cm -2 is reached, and the deeper 
region contributes only ~10% to the emergent emission. Correspondingly, the two O I flux 
ratios stay nearly constant between iVn = 10 22 and 10 25 cm -2 . On the other hand, the 
n(Ca ll)/n(0 I A8446) ratio becomes ~ 2 times larger when the c olumn density is increased 
from N-r = 10 23 to 10 25 cm" 2 . This effect was first pointed out by lFerland fe Persson! (119891 ) 
and is already taken into account in the above arguments. Hence, we concluded that the 
assumed column density does not significantly change the results from the standard model 
unless it is much smaller than ~10 22 cm -2 . Note that the O I line formation would be 
significantly suppressed in such a thin gas (see Fig. so that the EWs fall well below the 
observed values. 



3.3.2. Other Models: Microturbulence and Incident Continuum 



Some BELR models indicate a presence of the microturbulence wi thin the gas, which 
includes a large velocity gradient in a wind (e.g., iKonigl fe Kartid 1994T) and nondissipative 
magn etohydrodynamic waves in magnetically confined clouds (jReeslll987l ; lBottorff fe Ferland 
20001 ) . We performed another set of calculations with the microturbulence included, whose 
velocity is ftm-b- ft is known to have a large influence on the Fe II line strengths through a 
combin ed effect of the increased continuum fluo rescence rate and the decrease d line optical 
depth (INetzer fe Wills! Il983l ; IVerner et al.l Il999l ) . Using a strong dependency, IVerner et al. 
(120031 ) showed that t> tur b = 5-10 km s" 1 gives the mos t plausible values to a ccount for the 



observed Fe II and Mg II emissions. On the other hand, iBaldwin et al.l (120041 ) reported that 
the photoionized BELR models cannot reproduce both the observed shape and EW of the 
Fe II UV bump unless the considerable amount of microturbulent velocity corresponding to 
"Wturb > 100 km s _1 is included. Thus, we considered two models with the microturbulent 
velocity set to t> tur b = 10 km s _1 (model 2) and 100 km s _1 (model 3). The input parameters 
of these models are summarized in Table El 

The results in models 2 and 3 are rather similar to the standard model (Fig. [7]), while 
the whole patterns of contour are shifted a little toward a high-density regime as the micro- 
turbulent velocity is increased. It is thought to be due to the decreased optical depth in O I 
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A1304, which then reduces the electron population in the 3s 3 S° level and thus the collisional 
excitation rate of O I A8446. As a result, it has to be compensated by the increased particle 
density to reproduce the observed degree of collision rate. It is also noted that the predicted 
EWs of O I lines are moderately decreased on the whole (nn, U) parameter plane in model 
3. Since the difference in the wavelengths of Ly/3 (1025.72 A) and the pumped O I transition 
2p 3 P - 3d 3 D° (1025.77 A) corresponds to the velocity of 15 km s _1 , the Ly(3 fluorescence 
efficiency is reduced in such a turbulent gas. The best-fit U) parameters are determined 
in the same manner as in the standard model and listed in Table [3, which are U) = 
(10 1L5 cm" 3 , 10~ 2 - 5 ) in model 2 and (10 120 cm" 3 , 10" 2 5 ) in model 3. The predicted line 
strengths in these specific models are also listed in the table. 



Finally, we cha nged the incident continuum shape to be much harder as suggested by 
Korista et al.l (119971 ) , which has the form (T cut , a uv , a x , a ox ) = (10 6 K, —0.5, —1.0, —1.4). It 
was also adopted by Baldwin et al. (2004). The microtu rbulent velocity is set to t^urb = 

100 



km s" 1 , which makes this model identical to model 8 of Baldwin et al.l ( 120041 ) . Model 8 was 
considered to reproduce both the observed shape and EW of the Fe II UV bump observed in 
AGNs with strong Fe II emission, and was proven to be one of the successful models (see §4.41 
for more discussions). This is our model 4, whose input parameters are summarized in Table 
|6j Although the incident continuum shape is dramatically changed, the results are nearly 
identical to those in model 3, except that the predicted O I EWs are further decreased. The 
best-fit parameters are identical to model 3, (%, U) = (10 12,0 cm -3 , 10" 2 ' 5 ). They are listed 
in Table [7] along with the predicted line strengths. 



4. DISCUSSION 

4.1. Excitation Mechanisms of the O I Lines 

Our results indicate that the O I emission observed in quasars is formed through Ly/3 
fluorescence and collisional excitation in gas with a density of wh ~ 10 11 ' 5 cm -3 , illuminated 
by the ionizing radiation corresponding to the ionization parameter U ~ 10 -2 ' 5 . It can be 
tested by investigating the observed strengths of other O I lines, which should or should not 
be accompanied by the strongest three lines, A1304, A8446, and A11287, in each excitation 
mechanism. We list the observed flux ratios of these lines measured to date in Table [HJ 
Theoretical predictions of the ratios for the four possible excitation mechanisms, i.e., Ly/3 
fluorescence, collisional excitation, continuum fluorescence, and recombination, are collected 
from the literature and also listed. In the case of pure Ly/3 fluorescence, all the ratios listed 
in the table should be zero; actually, this is satisfied in most objects with one exception, 
1H 1934—063. The Ly/5 fluorescence plus collisional excitation from the 3s 3 S° level, which 
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our results indicate, would also predict all the ratios in the table to be nearly zero because 
the collisional excitation would predominantly produce the A8446 line photons against other 
lines listed. 



The detection of A7774 in 1H 1934-063 led iRodrieuez-Ardila et all (I2002ri ) to conclude 
that collisional excitation plays an important role in the O I line formation in NLSls. They 
stated that the observed A7774/A8446 ratio is ~0.2 after the removal of Ly/5 fluorescence 
contribut ion, whi c h is c onsistent with the theoretical prediction of pure collisional excitation 
given bv iGrandil Jl98oh (A7774/A8446 ~ 0.3; the "collisional excitation" row in Table [8J. 
However, the prediction was derived from the collisional cross sections only from the ground 
term (2p 3 P). It should be an oversimplification because the collisional excitation from the 
excited 3s 3 S° level, which we have shown is a large contributor to the A8446 emission, is 
not included. Then the observed strength of A7774 in 1H 1934—063, as large as 20% of 
A8446 (after Ly/3 fluorescence contribution was subtracted), should be attributed to other 
excitation mechanisms. Since another line, A7990, is also detected in similar strength in 1H 
1934—063, recombination migh t be working to some extent and producing these emissions, as 
Rodrfguez-Ardila et al.l (l2002bl ) recognized. An absence of the fourth mechanism, continuum 
fluorescence, is evident from the observed weakness of A13165; it is detected in none of 
the observed samples, while it should hay e a strength comparable to A11287 in the pure 
continuum fluorescence case (jGrandilll980D . 



Yet another mechanism of the O I line formation was suggested by IGrandil (119831 ). in 
which the 3s 3 S° level could decay to the metastable terms of the ground configuration via 
the semiforbidden lines A1641 and A2324 (see Fig. [5]). He suggested that fully half of the 
A 1304 photons wou ld be destroyed by this mechanism in the typical BELR cloud. However, 
Laor et al.l (jl997bl ) reported that they found no strong line at 1641 A in the spectrum of I 
Zw 1 and that such a line cannot add more than ~30% to the observed A1304 flux. 



4.2. Constitution of the A1304 Bump 



The A 1304 bump typically seen in AGN spectra is the result from a blending of an O I 
tri plet at A1302 . 17, A1304.86, and A1306.03, with a mean laboratory wavelength of 1303.50 
A JVerner et al.lll996l : IConstantin et alil2002h and a Si II doublet at A1304.37 and A1309.27 
(mean laboratory wavelength 1307.63 A, hereafter referred to as the A1308 doublet). This 
blending is very strong in most quasars, preventing a reliable estimate of relative contribu- 
tions of O I and Si II to the blend. Accordingly, deblending of the bump has be en tried mostly 



for NL Sls and quasars having relatively narrow broad emission line profiles. iBaldwin et al. 



( 119961 ) performed the method to fit the synthetic spectra constructed separately for each 



-16- 



quasar to deblend the bump. Although they could not firmly determine the strengths of 
both the O I and Si II components for any single quasar, the fitting results show that the O I 
fractions in the bump vary in the range f r om <3 0% to >60%. More plausible measurements 
were achieved for I Zw 1 by lLaor et all (Il997bl ). The Si II A1309.27 component is clearly 
detected in I Zw 1, which enables the reliable estimate of the O I and Si II fractions in 
the blend. Their template fit shows that 50% (56%) of the bump flux is due to O I in the 
optically thick (thin) Si II doublet case. iRodriguez-Ardila et al.l (j2002bl ) applied the same 
method to three NLSls and found that the average portion of the O I flux is 75%. 

Now we investigate the composition of the A1304 bump from the standpoint of the pho- 
toionized model calculations. We are particularly interested in whether the dominant con- 
tri butor to the bump, i.e ., either O I or Si II, actually varies from quasar to quasar as shown 
bv iBaldwin et all Jl996h . We plot the O I n(A1304)/n(A8446) versus n(A11287)/n(A8446) 
relation in Figure [9], which is calculated in models 1-4 with the (hh, U) parameters around 
the best-fit values with which the O I line formation is dominated by Ly/3 fluorescence and 
collisional excitation, as revealed above. The range of the observed n(A11287)/n(A8446) 
value is shown by two dotted lines. In order to obtain a rough estimate independent of the 
selected model, we regard the predicted n(A1304)/n(A8446) values enclosed by two dashed 
lines in the figure, drawn to include most of the points in the observed n(A11287)/n(A8446) 
range, as the allowed values against n(A11287)/n(A8446). These lines are expressed as 



n(A1304)/n(A8446) = 1.15 x n(A11287)/n(A8446) - (0.20 ±0.05). (1) 

Then the observed n(A11287)/n(A8446) ratios and the bump (O I + Si II A1304) flux give 
the estimates of the O I fraction in the bump; the results are listed in Table [9j It shows that 
the O I fractions actually vary from quasar to quasar by a significant amount, ranging from 
<20% in QSO J1139-1350 to >60% in 3C 273 and QSO B0850+440. While it was hoped 
that O I A1304/A8446 ratio would be a reliable reddening indicator of AGNs a s long as its 



intrinsic value was accurately determined (e.g., IRodriguez-Ardila et al.l l2002bl ) . our results 
show that it is more difficult than had been expected without accurate line deblending or 
precise model calculations. 

The derived fluxes of the Si II A1308 doublet are comparable or even a few times larger 
than those of O I A1304. However, it cannot be reproduced in our model calculations; the 
equivalent width of Si II A1308 is at most 3 A on the whole (nn, U) parameter plane in all 
models while EW (O I A1304) ~ 100 A (10 A if the covering factor is set to 0.1) with the 
best-fit parameters. Changing the assumed column density does not settle the problem, in 
which the Si II A 1308 flux is enhanced only by ~20% when N-r is increased from 10 23 to 
10 25 cm -2 . A possible explanation comes from the LOC scenario. We plot the distribution 
of EW (Si II A1308) as a function of the gas density rin and the ionization parameter U 
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in the standard model in Figure [TO], which shows that the Si II A1308 emission arises from 
the vast area on the parameter plane in similar strengths, while O I emission arises from 
the relatively limited region (Fig. [71 top right). Thus, the LOC-type integration over the 
parameter plane could significantly strengthen the predicted Si II emission relative to O I, 
which could account for the observations. An alternative way to interpret the strong Si II 
A 1308 emission is an unusually large microturbulent velocity present in the emitting gas, 



as suggested by iBottorff et all (120001 ) . According to their calculations, Si II A1308 is one 
of the most enhanced lines by the increased micoruturbulent velocity. They showed that 
very large values of the mi croturbulent yelocity are required in order to reproduce the Si II 



line strengths reported by Baldwin et all (119961 ). which correspond to f tur b > 1000 km s 1 



Si II A1308 is enhanced by more than 2 orders of magnitude in such a highly turbulent gas 
compared to the no microturbulence case, while O I A1304 is reduced by more than 80% due 
to the significant suppression in the Ly/3 fluorescence rate. 



Still other mechanisms of the Si II A 1308 enhancement are discussed in detail in Baldwin et al 



( 119961 . Appendix C). They suggested that dielectric recombination or charge transfer might 
produce significant emission in the highly excited Si II lines, including A1308, in order to 
account for their observed strengths relative to the lowest line, A1814, in one of the gas com- 
ponents (dubbed "component A") in the quasar Q0207— 398. Note that component A was 
found to have an extremely high (tih > 10 12 ' 5 cm -3 ) gas density, which they stated makes 
the above processes competitive through the thermalization of all resonance lines. 

Finally, it should be noted that the intrinsic values of O I + Si II A1304/O I A8446 
ratio could be larger than measured by more than an order of magnitude due to the internal 
reddening effect. Then the O I fractions in the A1304 bump would become even smaller 
than derived here, corresponding to the applied amounts of dereddening. An additional but 
probably smaller uncertainty comes from the possible variability event of the emission lines 
(see §£5p . 



4.3. Situation in Other Type 1 AGNs 



Observations of th e O I flux ratio between A8446 and A11287 in Syls and NLSls hay e 
been reported to date (ILaor et al.lll997bl ; iRudy et al.ll2000l ; iRodrfguez-Ardila et al.ll2002bl ). 
which are listed in Table [5] (see the Syl and NLS1 sections). The ratio of 0.55 ± 0.08 in the 
Syl is very similar to those of quasars, which may indicate similar emitting gas properties 
in the BELR of both types of AGNs. However, a larger sample of Syls is obviously needed 
in order to make a significant comparison. 
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On the other hand, the mean O I n(A11287)/n(A8446) ratio observed in NLSls is 0.73 
± 0.20, which is marginally higher than those of quasars and Syls (although they agree 
within the 1 a deviation). This situation is illustrated well in Figure [TTJ, in which we plot 
the number distributions of the ratio separately for quasars plus Syls and for NLSls. It is 
direct evidence that radiation-related processes are more efficient over collisional processes 
in the BELR gas in NLSls than in quasars and Syls. A possible explanation might come 
from the general idea of the NLS1 containing less massive black holes with higher mass 



accret i on rates, radiating at close to Eddington luminosity (IPounds et al.l 1 19951 ; lLaor et al. 



1997al ; iMineshige et al.ll2000i ). than normal Seyfert galaxies. Consequently, the line-emitting 
gas might be exposed to stronger continuum radiation so that the radiation-related processes 
become predominant. However, in the LOC scenario, the line-emitting region simply moves 
farther out when the continuum flux is increased at a given radius, so that the incident 
continuum flux (presumably the ionization parameter) in the region is kept n early constant. 



Such a situation is actually confirmed by the reverberation mapping results (IPeterson et al. 



20021 ). The difference in the ionizing continuum shape would also not be the origin of the 
O I n(A11287)/n(A8446) difference as we see in §3.3.2t in fact, the O I line ratios predicted 
in models 3 and 4 are quite similar to each other when the gas density and the ionization 
parameter are fixed (see Table Cj). 

Alternatively, the somewhat smaller gas density in the BELR of NLS1 gives a natural 
explanation for the decreased collisional excitation rate. In this case, the difference in the 
mean O I n(A11287)/n(A8446) ratio between quasars plus Syls and the NLSls corresponds 
to a few x 0.1 dex smaller gas density in NLSls, as seen in Figure [7] (top left; around the 
best-fit parameters log % ~ 11.5 and log U ~ —3.0). 



4.4. Role of the Relevant Gas in the BELR 

Our results suggest that a rather high (n# ~ 10 11 ' 5 cm -3 ) density gas is present in the 
BELR at the region illuminated by the continuum radiation corresponding to the ionization 
parameter U ~ 10~ 2 ' 5 , producing a dominant part of the observed O I and Ca II emissions. 
The gas is expected to be located at an outer portion of the BELR. It is a natural consequence 
of the low ionization potentials of the relevant ions relative to other typical emission lines 
observed in AGNs and is also indicated by the observed line widths. We list the full widths 
at half maximum (FWHMs) of the observed C IV A1549 and O I A11287 line profiles in 
Table [10] as the representatives of the high- and low-ionization lines. These two lines were 
chosen because they are relatively isolated from other lines and are strong enough to allow 
reliable measurements. The table shows that the line widths of O I A11287 are 0.3 - 0.5 times 
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narrower than those of C IV A1549. Assuming that the BELR kinematics is Keplerian so 
that Vfwhm oc r -1 / 2 , where Vfwhm is the emission - line width and r is the d i stance between 
the line-emitting gas and the central source (e.g., iPeterson fe Wandellll999l . 120001 ) . the O I 



emitting gas seems to be orbiting at a radius several times larger than that of C IV A1549. 
The reverberation mapping results for the low- ionization lines such as Mg II, Fe II, an d H/3 
( jDietrich fe Kollatschnylll995l ; IPeterson fe Wandellll999l ; IVestergaard fe Peterson! 120051 ) also 
support this picture, although the results for the relevant O I and Ca II lines have not been 
obtained to date. 

In order to check the consistency of our results with the observations of other emission 
lines, we list the calculated EWs of typical UV to optical lines seen in AGNs in Table [TTJ 
The calculation was performed with the best-fit (tt-h, U) parameters to the O I and Ca II 
observations and the covering factor of 0.1 in the standard model, while similar results 
were obtained in other models. We also list the EWs m easured in a Large Bright Quasar 
Survey (LBQS) composite spectrum (IFrancis et al.lll99ll ) for the compared observation. The 
composite spectrum was constructed from 688 quasars (absolute magnitudes Mb, < —21.5) 
and 30 AGNs (absolute magnitudes —21.5 < Mbj < —20.5) showing no strong evidence for 
the presence of broad absorption line troughs (see iFoltz et al.l 119871 ; iMorris et all Il99ll . for 
details of the observing procedures). The table shows that the O I and Ca II emitting gas 
exclusively contributes to the permitted low-ionization lines such as Mg II, UV Fe II, and 
H/3, as we expected. On the other hand, a main part of high-ionization lines such as Si IV 
+ O VI] A 1400 and C IV A 1549 is not produced in this gas; these lines should be formed 
in an inner portion of the BELR where gas is illuminated by stronger ionizing radiation. 
The forbidden line [O III] A5007 is also deficient in the emitted spectrum from the gas, 
which is consistent with the general picture that the line is formed in a much lower density 
environment in the NELR. 

The shape of the Fe II UV bu mp is correctly repro duced in models 3 and 4, in which 
fturb is set to 100 km s -1 as found by lBaldwin et al.l (12004 ). In fact, the best-fit parameters in 
these models, (%, U) = (10 12,0 cm -3 , 10~ 2,5 ), lie on the parameter region consistent with the 
observed EW and "spike/gap" ratio they defined (see their Fig. 7; our models correspond to 
log tt-h = 12.0 and log $h = 20.0). Models 1 (the standard model) and 2 cannot reproduce the 
observed spike/gap ratio; the bump is dominated by the strongest resonance lines, making 
the predicted ratio much larger than observed. On the other hand, none of the models can 
reproduce th e strong optical Fe I I emitters. It has long been a c hallenge to the photoi onized 



models (e.g.. ICollin fe Jolyl 120001 ) and was recently reviewed by Baldwin et al.l (120041 ) . who 



stated that the bulk of the optical Fe II emission might come from a region different from 
that which produces the UV Fe II e mission, such as collis ionally ionized gas. We will return 
to this issue in a forthcoming paper (ITsuzuki et al.l 120071 ) . 
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4.5. Interpretation of the results in the LOC Scenario 

In the context of the LOC scenario, an observed spectrum is the sum of emission from 
the gas at a range of different distances from the central source and also within a range of 
different densities. Hence the observed emission line flux should be expressed as 

-Fiine oc J J r 2 F(r,n u ) f{r) g(n R ) dn u dr, (2) 

where F(r, tie) is the emission-line flux of a single cloud at radius r from the central source 
and with density %, while /(r) and g(nu) are distribution functions of the gas covering 
fraction and density, respectively. 

In this scenario, the best-fit (%, U) parameters we derived above correspond to the 
model producing a dominant part of the observed O I and Ca II emission. We plot the 
quantity 

DF(r, n H ) = r 2 F(r, n H ) f(r) g{n n ) n u r, (3) 

for O I A8446 in Figure d2] to see which region on the parameter plane actually dominates 
the emission. The covering fraction and the gas-density distributions are assumed to be 
[f(r) oc r°, r -1 , r~ 2 ] and [g{n-a) oc n^, , %[ 2 ], respectively. With g(nn) oc or 
n^ 2 density distributions, the O I A8446 emission is dominated by much higher or lower 
density gas than the best-fit model. On the other hand, an f(r) oc r° type covering fraction 
results in dominant emission from the gas that is much farther from the central source (i.e., 
illuminated by much weaker ionizing radiation) than the best-fit model, while /(r) oc r™ 2 
predicts the opposite end. Thus, the functions of /(r) oc r _1 and g(nn) oc n# represent 
the most plausible distributions of the BELR gas. Note that it is also true in other models 
(models 2-4). It is a consequence of the fact that any of the assumed (%, U) parameters do 
not predict the EW of O I lines much greater than observed, i.e., EW (O I A8446) ^> 10 A, 
with a reasonable covering factor (see Fig. top right), so the best-fit parameters are chosen 
near the EW (O I A8446) maximum on the parameter plane. Since the emission-line flux of 
a single cloud can be written as F(r, tih) = (EW)F cont i nU um oc (EW)r -2 , with f{r) oc r~ l 
and oc n^ 1 distributions, DF becomes 

DF(r, n H ) oc r 2 [(EW)r~ 2 ] r" 1 n^ 1 n H ra EW. (4) 

It means that the DF distribution is identical to those of EW on the parameter plane so 
that the best-fit (%, U) parameters chosen near the EW maximum should also be near the 
DF maximum. On the other hand, other shapes of f(r) and g(nn) functions predict the 
DF maxima very differently from that of EW, i.e., very differently from the best-fit (nn, U) 
parameters. Thus, our results point to the gas distribution of f(r) oc r^ 1 and oc n^ 1 

forms in the LOC model to a first-order approximation. 
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On the other hand, the DF maximum for the Ca II lines occurs at a much lower value 
of U than the best-fit (tt-h, U) parameters with f{r) oc r _1 and g(nn) oc n^ 1 distribution 
functions, as shown in Figure US It implies that the gas distribution should have a cutoff 
somewhere around U ~ 10 -3 ' 5 at nn = 10 110 — 10 120 cm" 3 , which should originate from the 



physical nature of the BELR. In fact, iBaldwin et al.l (119951 ) constrained a range of the LOC 
integration over the parameter space to $ > 10 18 s" 1 cm" 2 , since gas at a larger distance than 
this limit w ill form graphite grai ns (jSanders et al.lll989l ) so that the line emission is heavily 
suppressed ( iNetzer &: Laorlll993l ). This limit corresponds to U > 10~ 35 at nn = 10 115 cm -3 , 
which is consistent with the above expectation. 



5. SUMMARY 

We present results of the near-IR spectroscopy of six quasars with redshifts ranging from 

0. 158 to 1.084. Combined with the UV data taken with the HST FOS, we give the relative 
strengths of O I A1304, O I A8446, O I A11287, and the Ca II triplet (A8498, A8542, and 
A8662). Detailed model calculations were performed in the framework of the photoionized 
BELR gas for a wide range of physical parameters and compared to the observations, which 
led us to the following conclusions. 

1. A rather dense gas with density nn ~ 10 11 ' 5 cm -3 is present at an outer portion of 
the BELR, whose distance from the central source corresponds to the ionization parameter 
U ~ 10~ 2 ' 5 . The gas is primarily responsible for the observed O I, Ca II, and Fe II lines 
based on the resemblance of their profiles. 

2. The strongest O I lines observed in typical AGN spectra, i.e., A1304, A8446, and A11287, 
are formed through Ly/3 fluorescence and collisional excitation. It is consistent with the 
absence or weakness of other O I lines such as A7774 and A13165. 

3. The fractions of the O I A1304 triplet and the Si II A1308 doublet blended in the A1304 
bump typically seen in AGN spectra vary from quasar to quasar, ranging from 20% to 
60% in the O I fraction. Such a large flux in the Si II doublet might be accounted for by 
the LOC-type integration over the gas physical parameter space or by an unusually large 
microturbulent velocity. 

4. The flux ratio between O I A8446 and A11287 observed in Syls may indicate that the 
physical condition within the line-emitting gas is similar to those in quasars. On the other 
hand, the ratios observed in NLSls suggest that radiation-related processes are more efficient 
over collisional processes than in quasars and Syls. It might be a consequence of a stronger 
continuum radiation incident on the line-emission region, or lower gas density in NLSls. 
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5. In the LOC scenario, the revealed physical parameters of the line-emitting gas should 
represent the maximum contribution to the observed O I and Ca II emission in the parameter 
space. Such a condition requires the distribution functions of the gas covering fraction and 
the density of f(r) oc r _1 and oc n^ 1 , respectively. 

We are grateful to an anonymous referee for giving a number of useful comments to 
improve this paper. G. J. Ferland kindly gave us the useful advices regarding the photoion- 
ization code Cloudy. We thank the staff of UKIRT for technical support and assistance 
with the observation. Use of the UKIRT is supported by the National Astronomical Ob- 
servatory of Japan. This work has been supported in part by Grants-in-Aid for Scientific 
Research (15253002, 17104002) and the Japan-Australia Research Cooperative Program from 
the Japan Society for the Promotion of Science. 
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Table 1: Characteristics of the Observed Quasars 



Object 


Redshift a 


M B b 


A v c 


J (2MASS) d H (2MASS) d 


J e 


H c 


3C 273 


0.158 


-26.9 


0.007 


11.8 


11.0 


11.5 




QSO B0850+440 


0.514 


-26.1 


0.011 


15.2 


14.4 




15.0 


3C 232 


0.530 


-26.7 


0.005 


14.9 


14.4 




14.4 


QSO J1139-1350 


0.560 


-26.3 


0.013 


15.5 


14.9 




15.0 


PG 1148+549 


0.978 


-27.7 


0.004 


14.5 


14.2 




14.6 


PG 1718+481 


1.084 


-29.8 


0.006 


13.5 


13.0 




13.0 



"Rcdshift measured in the UV spectra (jEvans fc Koratkarll2004l ) . 
& The B-band absolute magnitude taken from IVeron-Cettv &: Veronl (|2003[) . 
c Galactic extinction ( Schlegel et al. 1998 ); i?v = 3.08 was adopted (Pei 19921 ). 
d The J- and iJ-band magnitudes measured in the 2MASS observations. 
e The J- and iJ-band magnitudes measured in this work. 



Table 2: Observing Journal of the Near-IR Spectroscopy 



Target 


Exposure Time (s) 


Grism 


Air Mass 


Rationing Star 


3C 273 


2400 


IJ 


1.1 


BS 4708 (F8V) 


QSO B0850+440 


3840 


JH 


1.5 


BS 3451 (F7V) 


3C 232 


7680 


JH 


1.2 


BS 3625 (F9V) 


QSO J1139-1350 


6720 


JH 


1.3 


BS 4529 (F7V) 


PG 1148+549 


2880 


HK 


1.2 


BS 4761 (F7V) 


PG 1718+481 


1920 


HK 


1.5 


BS 6467 (F4V) 



Table 3: Observing Journal of the Near-IR Imaging 



Target 


Exposure Time (s) 


Filter 


Air Mass 


Standard Star 


3C 273 


50 


J 


1.1 


FS 132 


QSO B0850+440 


300 


H 


1.3 


FS 125 


3C 232 


300 


H 


1.6 


FS 127 


QSO J1139-1350 


300 


H 


1.6 


FS 129 


PG 1148+549 


300 


H 


1.2 


FS 121 


PG 1718+481 


150 


H 


1.6 


FS 141 



Table 4. Measured Emission-Line Fluxes a 



Emission Lines 


3C 273 


QSO B0850+440 


3C 232 


QSO J1139-1350 


PG 1148+549 


PG 1718+481 




UV Lines 


Ly/3 + O VI A1030 


2530 ± 140 


160 ± 24 




123 ± 10 


97.1 ± 9.9 


123 ± 9 


Lya + N V A1220 


19600 ± 200 


711 ± 26 


457 ± 16 


1030 ± 20 


1200 ± 20 


2070 ± 30 


Si II A1264 


53.5 ± 11.2 


5.88 ± 0.55 


3.60 ± 1.66 


<7.7 


14.4 ± 6.3 


25.6 ± 2.3 


O I + Si II A1304 


404 ± 50 


24.3 ± 3.0 


24.7 ± 3.0 


19.2 ± 2.0 


37.3 ± 4.2 


45.0 ± 8.2 


C II A1335 




10.2 ± 1.3 


9.94 ± 1.12 


7.01 ± 1.68 


24.1 ± 3.2 


20.3 ± 3.7 


Si IV + O VI] A1400 




70.1 ± 3.4 


41.1 ± 2.36 


68.6 ± 7.4 


101 ± 6 


174 ± 14 


C IV A1549 


7400 ± 140 


263 ± 7 


165 ± 6 


578 ± 7 


216 ± 8 


196 ± 31 


C III] A1909 


1680 ± 90 


112 ± 7 


128 ± 5 


51.9 ± 3.3 






Near-IR Lines 


O I A8446 


239 ± 26 


5.07 ± 1.17 


0.98 ± 0.17 


4.31 ± 0.77 


8.10 ± 0.26 


9.62 ± 1.05 


Ca II A8579 b 


343 ± 31 


8.15 ± 1.83 


<0.39 


2.08 ± 0.59 


9.40 ± 0.29 


10.5 ± 1.7 


Pae +[S III] A9540 


228 ± 12 


4.49 ± 0.54 


5.61 ± 0.20 


5.60 ± 0.47 


6.51 ± 0.35 




Pa<5 + Fc II A10050 


302 ± 7 


5.38 ± 0.37 


2.35 ± 0.22 


10.5 ± 0.8 


9.01 ± 0.29 


17.3 ± 0.48 


Fe II A10501 




1.95 ± 0.28 




1.25 ± 0.35 


0.99 ± 0.13 




He I A10830 


499 ± 20 


14.2 ± 1.4 


15.0 ± 1.3 


23.2 ± 3.6 


23.9 ± 13.1 


25.6 ± 1.9 


Pa7 A10941 


402 ± 20 


4.52 ± 1.15 


9.76 ± 1.27 


11.5 ± 3.6 


16.5 ± 13.1 


35.9 ± 2.1 


O I A11287 


59.4 ± 1.2 


2.34 ± 0.20 


1.29 ± 0.57 


0.77 ± 0.11 


1.99 ± 0.10 


3.41 ± 0.76 



a The fluxes are given in units of 10 15 erg s 1 cm 2 . 

b Total flux of the near-IR Call triplet (A8498, A8542, and A8662). 



Table 5. O I and Ca II Line Strengths Used to Constrain the Models 



Object EW (O I A8446) (A) O I n(A11287)/n(A8446) n(Ca II)/n(0 I A8446) a References' 3 



Quasar 



3C 273 


19.9 


0.33 ± 0.04 


1.46 ± 0.21 


This work 


QSO B0850+440 


37.7 


0.62 ± 0.15 


1.63 ± 0.53 


This work 


3C 232 c 


2.7 


1.76 ± 0.83 


< 0.40 


This work 


QSO J1139-1350 


25.0 


0.24 ± 0.05 


0.49 ± 0.16 


This work 


PG 1148+549 


51.5 


0.33 ± 0.02 


1.18 ± 0.05 


This work 


PG 1718+481 


15.9 


0.47 ± 0.12 


1.11 ± 0.22 


This work 


PG 1116+215 


23.8 


0.76 ± 0.27 




1 


Syl 


NGC 863 




0.55 ± 0.08 




2 


NLS1 


1H 1934-063 




0.64 ± 0.05 




2 


Ark 564 




0.82 ± 0.03 




2 


Mrk 335 




0.64 ± 0.05 




2 


Mrk 1044 




0.42 ± 0.05 




2 


Ton S180 




1.08 ± 0.16 




2 


I Zw 1 




0.76 ± 0.11 




3, 4 



a The representative wavelength of 8579 A was used to convert the Ca II triplet flux to the photon number flux. 

b References — (1) iMatsuoka et alj l l2005h . (2) iRodriguez -Ardila et al. | l|2002bh . (3) lLaor et all <1997Vjh . (4) 
iRudv et all fcooch . 

c This quasar is excluded from the sample to be compared with the model calculations, due to its anomalously 
weak O I A8446 emission. 
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Table 6: Input Parameters of the Photoionized Models 



Model 


Continuum 


Shape 




"^turb 




(T cut [K], a uv , 




) a 


(km s _1 ) 


1 (standard) 


(1.5 x 10 5 , -0.2, 


-1.8, - 


■1.4) 





2 


(1.5 x 10 5 , -0.2, 


-1.8, - 


■1.4) 


10 


3 


(1.5 x 10 5 , -0.2, 


-1.8, - 


■1.4) 


100 


4 


(1.0 x 10 6 , -0.5, 


-1.0, - 


■1.4) 


100 



a See text for the definition of these parameters. 



Tabl e 7: Results of the Photoionized Model Calculations and the Compared Observatio n 
Parameter/Prediction Model 1 Model 2 Model 3 Model 4 Observed 



Best-Fit Parameter 



n H (cm 3 ) 


10 11.5 


10 n.5 


10 i2.o 


10 120 


N/A 


U 


1Q -3.0 


1Q -2.5 


10 -2.5 


10 -2.5 


N/A 


Predicted Line Strengths 


EW (O I A8446) a 


12.0 


12.1 


8.72 


5.08 


> 10 


O I n(A11287)/n(A8446) 


0.56 


0.58 


0.52 


0.53 


0.46 ± 0.18 


O I n(A1304)/n(A8446) 


0.45 


0.48 


0.40 


0.41 




n(Ca II)/n(0 I A8446) 


2.51 


1.57 


0.95 


1.04 


0.1 - 5.0 



a The EWs are computed with a covering factor of 0.1 and given in units of angstroms. 



Table 8. Observed and Theoretical Values of the O I Line Ratios 



Object 


type 


A7254/A8446 


A7774/A8446 


A7990/A8446 


A13165/A11287 


References' 1 




Observed 


NGC 4151 


Syl 


< 0.08 


< 0.02 


< 0.04 




1 


NGC 5548 


Syl 


< 0.05 


< 0.06 


< 0.02 




1 


Mrk 42 


Syl 


< 0.03 


< 0.08 


< 0.03 




1 


Mrk 335 


Syl 


< 0.01 


< 0.04 


< 0.02 




1 


Akn 564 


Syl 


< 0.04 


< 0.03 


< 0.02 




1 


1H 1934-063 


NLS1 


< 0.03 


0.07 


0.04 


< 0.02 


2 


Ark 564 


NLS1 








< 0.07 


2 


Mrk 335 


NLS1 








< 0.05 


2 


Mrk 1044 


NLS1 








< 0.14 


2 


Ton S180 


NLS1 




< 0.08 


< 0.09 




2 


Theoretical Predictions 


Ly/3 fluorescence 




0.0 


0.0 


0.0 


0.0 




Collisional excitation 






0.3 b 






1 


Continuum fluorescence 




0.025 




0.052 


~ 1.0 


1 


Recombination 






7.2 C 






This work 



a References — (1) Grand 3 jl980h . (21 iRodn'guez-Ardila etaH <2002bh . 

b Only the collisional cross sections from the ground term are considered (see text). 

c Appropriate in the low-density limit (e.g., in the normal Galactic gaseous nebulae). 
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Table 9: Constitution of the A1304 Bump 3 



Object 


I b 


Si II 


3C 273 


264 (0.65) 


- 396 (0.98) 


8 - 


140 


QSO B0850+440 


15.3 (0.63) 


- 18.5 (0.76) 


5.8 


- 9.0 


QSO J1139-1350 


1.7 (0.09) 


- 4.1 (0.21) 


15.1 


- 17.5 


PG 1148+549 


7.8 (0.21) - 


12.6 (0.34) 


24.7 


- 29.5 


PG 1718+481 


19.8 (0.44) 


- 25.4 (0.56) 


19.6 


- 25.2 



"Fluxes are given in units of 10 15 ergs s 1 cm 2 , as in Table @] 
^Values in parentheses represent the fractions in the bump. 



Table 10: Measured FWHMs of the C IV A1549 and Q I A11287 Lines 



Object 


C IV A1549 


I A11287 


3C 273 


3940 


2160 


QSO B0850+440 


5540 


1960 


3C 232 


7430 


2410 


QSO J1139-1350 


3790 


1210 


PG 1148+549 


5860 


1830 


PG 1718+481 


3720 


1840 



Note — The FWHMs are given in units of kilometers per second. 
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Table 11: Equivalent Widths of the Typical UV and Optical Emission Lines 



Line 


Model a 


Observed 13 


Lja A1216 + N V A1240 


40.8 


52 


Si IV+O VI] A 1400 


0.49 


10 


C IV A1549 


1.35 


37 


C III] A1909 


0.23 


22 


Fe II AA2240 - 2650 


16.7 


49 


Mg II A2800 


65.0 


50 


H/3 A4863 


31.0 


58 


[0 III] A5007 


0.0020 


15 



Note — Equivalent widths are given in units of angstroms. 
a The standard model with the best-fit (njj, U) parameters to the O I and Ca II observations and a covering 
factor of 0.1. 

& The EWs of Emission lines in the LBQS c omposite spectrum we re taken from lFrancis et a~ (1991), except 
for that of the Fe II, which was taken from iBaldwin et ah ( 2004 ). 
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< 0.0 



3C 273 



Poe + [Sill] Hel A10830 




01 M 1287 



Fell A 1.0501- ■ 



9000 10000 11000 12000 

Wovelength (A) 

QSO B0850+440 



13000 



14000 




12000 14000 16000 

Wavelength (A) 

3C 232 



18000 



Poe + [Sill] Hel A10830 



Ol A1 1287 



12000 14000 16000 

Wavelength (A) 



18000 




Fig. 1. — Reduced near-IR spectra (solid lines) along with the sensitivity curves (with 
arbitrary amounts of scaling; dotted lines) for 3C 273, QSO B0850+440, and 3C 232 in the 
observed frame. Positions of the identified emission lines are marked with vertical lines. 
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QSO J1 139-1350 




12000 



o 0.5 b 
^ 0.0 



14000 16000 
Wavelength (A) 

PG 1 148+549 



18000 



Pae + [Sill] 



Hel A 10830 



Pay 

[ 01 A1 1287 




14000 16000 18000 20000 

Wavelength (A) 

PG 1718+481 



22000 



24000 



Pae + [Sill] Hel A10830 



01 A1 1287 




14000 16000 18000 20000 

Wavelength (A) 



22000 



24000 



Fig. 2.— Same as Fig. fflbut for QSO J1139-1350, PG 1148+549, and PG 1718+481. 
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Fig. 3.— UV spectra around O I A1304 for 3C 273, QSO B0850+440, 3C 232, QSO 
J1139-1350, PG 1148+549, and PG 1718+481 in the observed frame. Positions of the 
identified emission lines are marked with vertical lines. 
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Fig. 4.— Fitting results for O I A1304, O I A8446 + the Ca II triplet (A8498, A8542, and 
A8662), and O I A 11287 in 3C 273. Observed spectra are shown by thin solid lines while the 
determined continuum and the fitted Voigt profiles are shown by thick solid lines. The sum 
of these fitted components are also shown by a thick solid line for O I A8446 + the Ca II 
triplet. The wavelength ranges used in the fitting are indicated by vertical dashed lines. 
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Fig. 5. — Partial Grotrian diagram of an O I atom. Solid lines represent permitted tran- 
sitions, while dashed lines represent semi-forbidden transitions. Transition paths of the 
Ly/3-pumped electrons (corresponding to A1026, A11287, A8446, and A1304) are shown by 
thick solid lines. The six energy levels considered in the model calculations are shown by 
thick horizontal lines. 
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Fig. 6.— O I photon flux ratios n(A11287)/n(A8446) (top) and n(A1304)/n(A8446) (bottom) 
as a function of the gas density rin and the ionization parameter U, calculated in the standard 
model. The contours are linearly spaced. Arrows indicate the direction toward which the 
incident-ionizing continuum flux <£> increases. 
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Fig. 7.— O I n(A11287)/n(A8446) (top left), EW (O I A8446) in units of angstroms (top 
right), n(Ca II)/n(0 I A8446) (bottom left), and their composite plot (bottom right) as a func- 
tion of the gas density nu and the ionization parameter U, calculated in the standard model. 
In the bottom right panel, the three constraints from the observed O I n(A11287)/n(A8446), 
EW (O I A8446), and n(Ca II)/n(0 I A8446) values are shown by thick solid lines, thin solid 
lines, and the shaded area, respectively The models consistent with all the constraints are 
marked with diamonds. Dashed lines and dotted lines represent mean values of the observed 
O I n(A11287)/n(A8446) and n(Ca ll)/n(0 I A8446), respectively. Note that the contours in 
the top left panel are linearly spaced, while they are spaced logarithmically in the top right 
and the bottom left panels. 
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Fig. 8. — Cumulative flux of O I A8446 (with an arbitrary scaling; solid line) and ratios of O I 
n(A1304)/n(A8446), I n(A11287)/n(A8446), and n(Ca ll)/n(0 I A8446) (dashed line, dot- 
dashed line, and dotted line, respectively), plotted against the gas column density measured 
from the illuminated surface. The model parameters other than the column density N-r used 
in the calculation are identical to those of the standard model with (n H , U) = (10 115 cm -3 , 
10" 3 - ). 
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Fig. 9.— O I photon flux ratios n(A11287)/n(A8446) vs. n(A1304)/n(A8446) calculated in 
models 1-4 with the (%, U) parameters around the best-fit values, with which the O I line 
formation is dominated by Ljj3 fluorescence and collisional excitation (see text). Diamonds 
represent model 1 with different (nu, U) parameter values, while squares represent model 
2, circles represent model 3, and triangles represent model 4. Dotted lines show a range of 
the observed values, n(A11287)/n(A8446) = 0.46 ± 0.18. We regard the n(A1304)/n(A8446) 
values enclosed by two dashed lines, drawn to include most of the points in the observed 
n(A11287)/n(A8446) range, as allowed values given the n(A11287)/n(A8446) ratio. 
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Fig. 10. — EWs of the Si II A1308 doublet in units of angstroms, as a function of the gas 
density % and the ionization parameter U calculated in the standard model. A diamond 
represents the best-fit parameters to the O I and Ca II observations. 
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Fig. 11. — Number distributions of the O I n(A11287)/n(A8446) ratios observed in AGNs. 
The top panel shows the distribution for quasars and Syl, while the bottom panel shows the 
same for NLSls. 
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Fig. 12. — Distributions of the DF (see text for its definition) values for O I A8446 as a 
function of the gas density nn and the ionization parameter U calculated in the standard 
model. Assumed distribution functions of the gas covering fraction and density are indicated 
at the top of each panel. Diamonds represent the best-fit parameters to the O I and Ca II 
observation, while triangles represent the parameters predicting the DF maxima. 
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Fig. 13. — Distribution of the DF (see text for its definition) values for Ca II as a function 
of the gas density nn and the ionization parameter U calculated in the standard model. 
Assumed distribution functions of the gas covering fraction and density are f(r) oc r _1 and 
g(nu) oc n^ 1 . A diamond represents the best-fit parameters to the O I and Ca II observation, 
while a triangle represents the parameters predicting the DF maximum. 



